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The form of the resulting Feynnian propagators in a proposed local and gauge invariant QCD for 

massive fermions suggests the existence of indefinite metric associated to quark states, a property 

that might relate it with the known Lee- Wick theories. Thus, the nature of the asymptotic free quark 

states in the theory is investigated here by quantizing the quadratic part of the quark action. As 

opposite to the case in the standard QCD, the free theory does not show Hamiltonian constraints. 

The propagation modes include a family of massless waves and a complementary set of massive 

oscillations. The theory can be quantized in a way that the massive modes show positive metric 

and the massless ones exhibit negative norms. It is remarked that, since QCD is expected to not 

fT^ ' exhibit gluon or quark asymptotic states, the presence of negative metric massless modes does not 

constitute a definite drawback of the theory. In addition, the fact that the positive metric quark 

("^ ' states are massive, seems to be a good feature of the model, being consistent with the approximate 

^SJ ' existence of asymptotically free massive states in high energy processes. 

Q_|. PACS numbers: 12.38.Aw;12.38.Bx;12.38.Cy;14.65.Ha 



< 



o 



cd 



I. INTRODUCTION 



In a recent work [l| it had been motivated an alternative for the description of quarks interacting through gluons, 
which possible equivalence with massless QCD was also argued. The outcome was a theory with an action given by 
Oh' the massless QCD one, which also incorporates additional terms, one for each quark flavour. The theory constitutes 
^ \ a local and gauge invariant completion of the one discussed in Ref. 043, and, up to our knowledge, it also furnishes 

[^ I . a new local formulation for strongly interacting massive quarks. The new Lagrangian terms appearing could also 
be interpreted as possible counterterms, within a special renormalization scheme for massless QCD. That is, it was 

'"^ , pointed out that such terms could be present in an effective action for massless QCD. These terms do not eliminate 
^ ' the power counting renormalizability of the model, because the new quadratic part of the action modifies the behavior 

[if' , of the free quark propagator, which now gets a 1/p^ behavior at large momenta. The expansion is expected to be 

_j. ' employed in further studies directed to explore physical questions as: the possibility of a flavour symmetry breaking 

f— ^ I effect determining a hierarchical quark mass spectrum, and the properties of quark-quark interaction potential in 
the scheme. However, the new diagrammatic expansion, shows a conceptual feature which is convenient to be 
studied before further employing it in the above cited studies. This special property is that the new more convergent 

X^ , free propagator is given by the difference between a free massless quark propagator and a massive one. This 
structure is very similar to the typical ones appearing in the so called Lee- Wick theories, in which usual bosons 
propag ators are modified by substracting massive propagators in order to make the expansion more convergent 
[10l4l2l |. Such substractions, for the bosonic fields, although making the propagators decreasing more rapidly at 

^^ I large momenta, introduce states of negative norm in the space of states, which mass shells are defined by the poles 
5h ■ of the new propagators [13l Il4| . Thus, the quark modified propagator appearing in reference [l|, which is the 
difference between a massless and a massive one as described before, might in principle lead to the appearance 
of an indefinite metric in the quark sector of the state space of the theory. As remarked in Ref. [l3| , the presence 
of an indefinite metric could not directly imply a limitation in the physical applicability of the theory. In our case 
this possibility might be indicated by the fact that for QCD, it is conceived that there are not asymptotic states for 
the elementary quarks and gluons. Also, the negative metric states in the Lee- Wick theories can be compatible with 
an unitary scattering matrix by the appearance of decaying rates for these negative norm states due to radiative 
corrections. Thus, the presence of those modes on the proposed variation of QCD might not automatically define its 
lack of applicability. However, the existence of such modes, could occasionally lead to drastic modifications of the 
predictions at very high energy processes, where the perturbative expansion is expected to approximately work due 
to asymptotic freedom. Thus, to know the type of negative metric states which can appear in the proposed theory 
becomes an important question. It can be also noted that the gluon free action is already including an indefinite 
metric associated to the unphysical gluon components. However, since the gluon and quark sectors do not couple in 
the free part of the action, the gluon free theory quantization is identical as in the usual QCD. 

Therefore, in this work, we investigate the quantization of the free part of the quark action of the proposed theory. 
For this purpose we first consider a reparametrization of the classical action obtained in [1| in order to express the 
relativistic quark propagator in a form in which the massive states appears as related with poles associated to positive 
metric states and the massless ones on the contrary, as linked with negative metric ones |16{ . Next, the description of 



the classical propagation modes of the free Lagrangian equations is considered. Two sets of quark modes arise. One 
of them is identical to the one corresponding to massless quarks in QCD, and a second one is a set of massive states of 
polarizations being equivalent to the usual set of modes satisfying the Dirac's equation. The classical Hamiltonian 
and the canonical equations are also written. A helpful form of the Hamiltonian, to be further employed for the 
quantization process is obtained. It was used the specific form of the generalized Hamiltonian procedure presented in 
Ref. 15]. The consistency of the Hamiltonian and Lagrangian equations were checked. Further the quantization is 
considered by defining field operators as superpositions of the previously determined wave modes. The Hamiltonian is 
then expressed in terms of the creation and annihilation operators defining the massive and massless quark fields. The 
result, as usual was written as a linear combination of the number operators for each quark mode with the coefficient 
defined by the energy for the mode. This conclusion followed after choosing the form of the anticommutation relation 
among all the creation and annihilation operators of the fermion modes, that define the quantum states as showing 
positive norm for the massive quarks and negative ones for the massless ones. Central in determining the chosen 
selection of the creation and annihilation operators for the single particle states and the metric of the state vector 
space, was, on one side, the requirement of defining a bounded from below value of the energy of the system. The 
process to implementing these properties was also basic in defining the positive or negative metric character of 
the various quantized modes. The satisfaction of the same classical field equations by the field operators and their 
momenta at all times was shown. The quantum commutators between all the fields and their momenta satisfied the 
usual quantization rule of being given by the imaginary unit times the results of the classical Poisson brackets. 

The presentation proceeds as follows. In Section II, the general form of the action for the theory is presented and 
the notations employed are defined. Afterwards, the general action and the quark propagators are rewritten, to fully 
evidence the similarity of the propagator structure with the ones appearing in the Lee- Wick theories. This is done in 
a way suggesting the positive metric of massive states and negative metric for massless states. Section HI, presents 
the determination of the quark modes associated to the poles of the propagator. The Section IV considers the writing 
of the Hamiltonian equations for the Grassman quark fields and their canonical momenta. The equivalence of the 
canonical and Lagrangian equations is verified and a formula for the Hamiltonian in terms of the fields is derived. 
Finally, in Section V the fields are consistently quantized by properly choosing anticommutation rules for the creation 
and annihilation operators for all the modes, and selecting the positive metric for the massive quarks and a negative 
metric for the massless states. 



II. MODIFIED QCD AND THE LEE- WICK THEORIES 

The classical action defining the model introduced in Ref. [l| had the explicit form 



(1) 



in which the conventions for the various quantities coincide with the ones employed in Ref. [l7[. In details, the 
appearing fields are defined as follows 
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where / = 1, ..., 6 indicates the fiavour index. The expressions for the Dirac conjugate spinors and covariant derivatives 
are 
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in which the Dirac's matrices, SU{S) generators and the metric tensor are defined in the conventions of Ref. [17|, as 

{7^ i"} = 2g''^ [Ta, n] = I r'^n. 7° = /3, 1' = /? «^ j = 1, 2, 3, 
/ 1 00 
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Other useful definitions and relations for the coordinates are 

tt- — i/j - — - yjj . ttj I - — - itt- m iij . j-i .iX' /. ''-'Li ' — " ^ LLi/ ''-' 1 "-^ — ^ ■ 

The complex conjugate operation is symbolized by the superscript *, by also assuming that the conjugate operation 
of a quantity also implies the Hermitian conjugation of any operator entering its definition. For example, in the 
conjugate transposed fermion field ^l {x) = {^h{x)Y'* the operator structure entering the definition of '^h{x) should 
be substituted by its Hermitian transposed structure. In what follows the color and spinor indices of the fields and 
other magnitudes will be omitted in order to simplify the writing. The time variable t will be defined by t = x^ . 
The system will be assumed to be enclosed in a large cubic spatial box of volume V ^ on which periodic boundary 
conditions are imposed for all the fields. This define the spatial momenta as taken the usual discrete values. 

The quadratic in the quark field of flavour /, free action of the theory is given by 



5o,/ = dx ^f{x){i^>'d^ - X d')-^iix) 

= - fdx^fix)Af{d)^,{x) (8) 



in which >f is a constant with dimension of length, which the discussion in former works motivating the proposal of 
the action ([T|) in Ref. [ll|, suggested to be related with quark condensation effects. Those analysis were associated 
to the consideration of a fermion squeezed state as the vacuum leading to modified non local perturbative Wick 
expansion for massless QCD . They strongly suggested the local action form ([T]) as possible effective action of the 
massless QCD g-Q. 

The Feynman propagator 5*/ for the considered quark of flavour /, is given by the above deflned inverse of the 
kernel A/ and satisfles 

-{ij^'d^~>cd^)Sf{x-y)^5{x-y), (9) 

or in terms of its Fourier transform Sf{x) — J dp Sf{p) exp(— ip.x)), 

-ij^p^ + xp^)Sfip) = I. (10) 

In order to simplify the notation, in what follows we will use the same symbol / (or simply the unity 1) for the 
Kronecker delta in the color and spinor indexes, for both cases: the two component and the four component spinors, 
when no confusion can arise. 

Thus, the quark propagator for the quark of flavour / has the form 
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The free quark propagator shows poles of mass ?ti/ = — , which are inversely proportional to the parameter xr. It 
also present poles of zero mass. This propagator has the structure given in the last line of pT|) . which expresses it, in 
the here employed notation, as the usual massless quark propagator minus the Dirac propagator for a massive quark 
in which the sign of the mass had been reversed. As mentioned in the introduction, this structure closely resembles 
the Pauli Villars kind of regularization employed in the Lee- Wick theories, which might lead to the appearance of 
negative metric states in the resulting quantum field theories. The massless component of the propagator shows the 



standard form for massless QCD with positive metric, and the massive component has a changed sign for the usual 
Dirac component of the action (as evidenced by the large momentum limit of this term). This suggests that the 
quantization of the Lagrangian with the employed field parametrization, will lead to negative metric for the massive 
fermion sector of the asymptotic states. However, since the theory shows space time inversion invariance, the sign 
of the time which could effectively correspond with the physical time is not well defined at the current stage of the 
analysis of the theory. Then, this freedom can be considered as factor to employed in the process of constructing 
the theory. After adopting this point of view, it seems possible to redefine the fields in the action by performing the 
following change of variables and coordinates within the functional integral defining the Green's functions generating 
functional 



■k , ,, -—/k , 



*(a;) = *)(-a;') = '^fix'), "^{x) = *}(-x') = Wf{x'), 

A'^'ix) = A'^'i-x') ^ -A'^^^ix'), 

cix) = c(-x') = c'ix'), c{x) = c(-x') = c'ix'). 

After performing this changes in the general expression for the action ([1]) , its expression in terms of the new primed 
fields and coordinates, for afterwards returning to the previous notation for the coordinates and fields by removing 
the primes in all the magnitudes, it follows for S 



S = 
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The new expression for the free quark action takes the form 



Soj ^ dx */(x)(-i7''a^ - X d^)^f{x) = - j dx-^f {x)kf{d)^,{x), (13) 

which defines for the new propagator kernel 

Sf{p) - 
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where now the massive term has the form of the typical Feynman quark propagator for massive quarks and the 
massless one has a different sign, which suggests its association to negative metric states. In order to define the precise 
way in which the theory implements the introduction of a space of states with indefinite metric for quarks, in next 
sections we perform the explicit quantization of this new form of the action. 

III. CLASSICAL QUARK PROPAGATION MODES 

Let us determine in this section the quark propagation modes defined by the free fermion component of the action. 
These u{x) modes satisfies the Lagrange equations defined by ([13]), which are 

{-ij'^d^, - K d^) u{x) = 0, u{x) = I dq^u{q) exp(-g,,a;^). (15) 

They imply for the Fourier components u{q), the matrix equation 

i-^q^ + K q^) u{q) = 0. (16) 

This expression clearly shows that there exist massless as well as massive propagation modes. Let us separately 
consider in what follows the massive and massless solutions. 



A. Massive solutions 

The Fourier four momentum q associated to any of the waves being determined, will be expressed in the form 
q ~ (g", g*). Since we are searching for the massive modes, in the Lorentz frame in which the momentum vanishes, 
the momentum can be written in the form q = ((jr°,0), and thus, the equation (jl6|) reduces to 

qoi-J +xqo)u{q)^Q^qo[^ (>. go + 1)/ j U^'V ' ^ ^ 

which clearly allows to easily define in this frame the solutions for the separate cases in which the energy component 
go takes positive or negative values. Both cases are considered in what follows in this subsection. 

1. Positive energy waves 
Below, it will be convenient to define a mass parameter mf as 

K=—. (18) 

TO/ 

For this case of positive energy waves, in the selected rest frame, the energy qo, takes the positive value qo — nif and 
the positive energy spinor polarization gets the simple form 



^'^(;)'^'^(;)- 



Then, the massive positive energy polarization of the considered theory coincide with the ones associated to the 
positive energy modes of the usual solutions of the Dirac's equation for a massive particle. Therefore, the spinor 
polarizations in any Lorentz frame might be written in the form 



-(^) - J '-fl^r i ^-^ 0r],r^l,2, (19) 



~^.f = a'p\ e™^ {-f) = Jm'j +f^,q^ (e„^. {'f),f), 



in which the energies of the modes are given by £mf{p) = A/^'f + P^- The polarizations had been chosen to satisfy 
the normalization conditions 

Note that the normalization for the spinors is fixed in a slightly different way as the usual one. This was done in 
order to assure the normalization to one of the spinors in the large quantization volume for the usual spinor scalar 
product. 

2. Negative energy waves 

In a similar way, for the negative energy case qo — —mf the solutions of the equation (J17p leads to the spinor 
polarization of the form 



•<^'^\/^^grHr"V-'^' 



"^.^ = crV, e™/(^) = \/mj+t^,q= (-e™^(^),^). 



in which again the energy parameter e^^ (p) — \ rni + ~p ^ and the polarizations were chosen to satisfy the normal- 
ization conditions 

em^ ( P ) 

The four momenta of the negative energy waves is given by q = (— Emj (~P )j ~P )• Note that as usual, in order to assure 
orthonormality for different values of the spatial momentum p of the modes, the spatial momentum index in the 
Fourier expansion of the original wave u(x), had been changed in sign for the negative energy solutions. 

B. Massless solutions 

These solutions of the general equation ([TO)) , since having q^ — (g")^ — q^ = , satisfy the massless Dirac 
equation 



i 7^*9^ u{x) = 0, u{x) = dq u(g) cxp(-g^a;^), (21) 

which for the Fourier components u{q) take the form 

rq^ uiq) = 0. (22) 

In terms of the null four momentum q — (g^, q*) it can be written as follows 

(go7° - iWHq) = /3(<Zo - c.''q^)u{q) 

qol -~a .~^ \ ( Pi 



It is clear that by considering the two component spinors /3i and P2 as eigenfunctions of the helicity operator 17 .q , 
the positive and negative energy polarizations can be determined. Below we examine the two cases. 

1. Positive energy massless waves 
After defining the (3'^^{p), (i^^CY) as the eigenfunctions of o -P , given by [1^ 



^{~f) = p^ = (ni, 712,713), 



n I 



W\ 
the positive energy solutions of momentum q = (eo(~p ), p ) can be written in the form 



M) = ^-(f'^^-^^],l-+l,-l, (24) 
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The normalization of these modes satisfy 
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2. Negative energy massless waves 
The negative energy modes, defined hy q^ — —\q\ , have the polarization 



^(^) = 7^ = (ni,n2,n3), 

where, as it is convenient for assuring orthonormality of the modes for different spatial momentum index, the spatial 
Fourier momentum q was fixed to be equal to minus —p ■ The appearing two component spinors have the same 
already defined explicit expressions 

and the four Fourier momentum of the mode has the form q = (— eo(~p ), p) with eo(~P ) = |~P I- These polarizations 
obey the orthonormality relations 

VoHf) viit) - 5"', v',{-^) vl{-^) = 0. (28) 

IV. FREE QUARK FIELD QUANTIZATION 

In this section we firstly present the application of the generalized Hamiltonian formalism to the free part of the 
quark action of the modified QCD in the form (fT3|). This discussion will also permit to obtain a formula for the 
Hamiltonian of the system which is further employed to quantize the theory. 

A. Generalized canonical procedure 

After performing few integrations by parts, the classical action S in (1131) can be expressed as a functional of the 
quark fields and their first time derivatives in a form which determines for the Lagrangian L of the system 

L= I dx^m\x){-'-a^^ ^ + ^ Cdt^t-^ :^)P)^!{x), (29) 

V^ = ^M - "d^.a^ = (/,ai,a2,a3),vl/t(a;) ^ ^>*T(^^)^ (30) 

where the right or left arrows appearing over the derivatives indicate whether they are acting on the right or the 
left expressions respectively. It can be helpful to recall that when operators will appear, the complex conjugation 
operation, designed by the superindex "*" also includes the Hermitian transpose of the operator structure. 

The momenta associated to the, up to now, classical fields * and ^^ are defined in the generalized mechanics with 
anticommuting variables (See Ref. [l5[) in the form 

V i 

in which the "right" Grassman functional derivative employed is defined for functionals which depend of the field 
and their derivatives as functions of the spatial variables "af at a time slice defined by a instant t. This instant is 
given by the fixed time slice employed to perform the spatial volume integration defining the Lagrangian in (|29p . 
It is interesting that this fermion theory is regular, since it does not have Hamiltonian constraints. Then, all the 
velocities can be explicitly expressed in terms of the coordinates (fields) and their momenta as follows 

dt^\x)P = -CT:^\x)+n^{x)), (33) 

XT Z 

Pdt^ix) = --(n^t {x) + '-^{x)). (34) 

>c 2 
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With these expressions, the Hamiltonian at the defined time shce can be expressed as a functional of the fields and 
momenta in the way 



H = I dx^{n^{x)dt'^{x)+n^i{x)dt'^\x)) - L 
= I dx^{-{-n^{x) I3n^i{x) + --^\x) (3-^{x)- 

J K 4 



-^\x) p n^,{x) - -n^{x) p ■^{x))+ 

^-^\x)a''l^,^{x) + X *t(a;)^.^/3^(a;)}. (35) 

The expression for the Hamiltonian, allows to write the canonical equations of motion in the form jl5| 

dt^\x) = j^^-^-,H = l(n*(x) + l-^\x))p 
dtii^ix) ^-H-^ = ^n^{x)p-^^\x)p 

o^(x) 2k Ak 

+ m!\x)a''^, + KV^^\x)li, 

d,uM-) = -i/^ = -^p n,t(.) + ^/? *(.) 

+ i a*5,*(x) - K^V'^^ix), (36) 

which after eliminating the momenta in terms of the fields and their time derivatives lead to the Lagrange equations 

{-i -y^df, - K 9^)*(.t) = 0, 
in consistency with the starting ones in (jlSp . 

B. Quantization 

Let us consider in this section the quantization of the quark fields \I' and ^^. For this purpose it seems convenient 
to find a formula for the Hamiltonian H in terms of the fields and their time derivatives. In the case of the Dirac's 
equation, and thanks to its constrained canonical structure, the Hamiltonian can be expressed in terms of the fields 
^ and ^^ and their spatial derivatives. However, the here considered mechanical system is regular and the formula 
for the enenrgy should also contain the time derivatives of the fields. The expression for H can be rewritten after 
eliminating the momenta in psp by employing the canonical equations ([36]) 



ii= f dx^x dt'i'^x) (3 dt-i'ix)+i^''{x)a'd^^{x) + x'$^{x) ^.'^/3*(a;)} 
= I dx^-^\x){>c'^ t fi~^ t + i a'd, + K^^ ^ (3)-^{x) 

dx^^^{x)h{d)^{x), (37) 



in which, as it was remarked before, the time derivatives of the fields enter the definition of the kernel h. 

Let us now start the quantization procedure by considering the quark fields as superpositions of all the oscillator 
modes with coefficients defined by operators. They will be identified as creation and annihilation operators for all 
the modes. For this purpose it will be useful to separate the field as a sum of two components: one associate to the 
massless modes and another defined by the massive ones. The kind of the mode will be represented by the greek 



index k taking two values: or m. Then the total field will be written as follows 






c|(yK(^)exp(z(eo(^)t + y.-^))) 
Mo(a;) +«o(2:) (39) 



))- 



*m(a;) - ^51 H (&r(^)u™,(^)expH (e„,^.(^)t - ^.^ 

dl {tVm,, it) exp(* (e„, {-f)t + y .1^)) 
= M,„(2:) +i;m(x), (40) 

where the polarization vectors for all the massless or massive modes were before determined. The assignation of the cre- 
ation or annihilation nature for each of the eight operators (a;(p ), al{p),ci{p ) , c] (p) ,br{p) , blCp) , dr{~p) , dlCp)) 
is not yet defined at this level. It will be fixed afterwards, by the requirements to be imposed on the metric of the 
state space to define a bounded from below energy for the many particle states. Thus, up to now the symbols f over 
any of these quantities will only mean that it is the Hermitian conjugate of the operator, not that it is a creation one. 

One important point to note is that the expression p7p is defined by the spatial integral in a given fixed time 
slice. Thus the time derivatives of the fields entering, will be evaluated by using the time dependences of the classical 
modes. Therefore, for the sake of consistency it should be required that after quantization of this free theory, the 
expression of the Hamiltonian H should reproduce the classical time dependency of the fields through the unitary 
transformation defining the operators in the Heisenberg representation '^{Xjt) = exp(iiJ t)^(~af ,0) exp(— iTf t). 

In addition, the field will be also decomposed in the sum of its the terms associated to the positive energy modes 
u^{x), plus the terms defined by the negative energy modes v^{x), for each of the two values of the mass k, as follows 

^ix)= J2 "«(^)+ E "-(^)- (41) 

K.— 0,771 K.— 0,m 

In terms of this expansion the Haniihonian has the expression 

- E E { I dx^ui{x)h{dW>{x)+ I dx\i{x)h{d)v^>{x)+ 

dx^ui{x)h{d)v^,{x)+ f dx^vl{x)h{d)u^,{x)). (42) 



The procedures for the calculation of the appearing matrix elements of the kernel h{d) are illustrated in the 
Appendix A. It follows that they vanish when both entering modes have different u or w types and also when 
they have different masses. The calculation of the few non vanishing terms leads to the following formula for the 
Hamiltonian operator 



H = dx^¥{x) h{d) *(a;) 

= / dx^ul,^{x) h{d) Um{x) + / dx^vl^{x) h{d) Vm{x)+ 
(ix^uj(x) h{d) uq{x) + I dx^vl{x) h{d) vq{x) 



= Y.Y. (-^o(^) a\{-f)ai(-^) + e,(-^) cK^)4(^)) + 
^ i=±i 

E E (^™/(^) ^l(~^)^r{t) + e,nA~^) dl{-f)dr{t))- (43) 

■f r=l,2 
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This expression is central for the purpose of deciding about the assignation of the creation and annihilation 
qualities to the introduced operators and the type of indefinite metric of the Fock space of states. This formula 
indicates that for massive modes quantities, the usual assignation br{p) and dr{p) to annihilation and bl{p) and 
(ij(~p ) to creation operators, assures that the contribution of their created particles (over a vacuum state annihilated 
by all the 6r(~P ) and dr{p)) to the energy will be bounded from below. Then, let us adopt the following commutation 
relations for the massive type operators 

{bU-fi),hAt2)} = Srn Sr^r'S^''Htl,f2), (44) 

{dli-f,),dr,if,)} = Srn <5„'<5 ^ (^i , ^ 2) (45) 

{A,B} ^ AB + BA, 

where the number Sm = ±1 yet allows to select positive or negative metric quantization of the associated states 
for massive fermions [ij, [IJ. The function (J^-^'Cp i, p 2) represent the Kronecker Delta function, equal to one for 
equal arguments and vanishing when they differ. In what follows, the definition of the Dirac's Delta function of the 
coordinates will be employed which is related with the discrete momenta in the quantization box as follows 

S(D)(j^^ _ ^2) = 5] ^ exp(zy.(^i - 1^2)) 

= J^eM^t-{'^i~^2)). (46) 

For the case of the massless modes quantities, it is clear that in order to the energy be positive, given the fermion 
character of the fields, the assignation can be to consider a creation type for the operators a;(~p ), ci{p) and an 
annihilation nature for the operators aj{p), cj(~p ) . With this interpretation, and since the particles are fermions we 
will define the adopted here creation and annihilation operators ai{p),ci{p) and their Hermitian conjugates as 

c,i-f) = cl{f), dlit)=ciif). (47) 

assuming to obey the anticommutation relations 

{aji-fi),ai,i-f2)}^soSij,S'^''H-fi,lt2),{cl{'fi),ci,i-f2)} = soSu,S^''H-fi,f2), (48) 

where in addition, any anticommutator between any two elements pertaining to the set formed by the operators and 
their Hermitian transposes, vanishes, if the two operators pertain to different classes of the four types "a ","c ","& " 
and "d ". The number sq = ±1 in this case, also allows for fixing positive or negative norm for the massless single 
particle states. Therefore, in what follows, the quantities ai{jJ), ci{~p), hr(p) and dr{j}) will be assumed to define 
the annihilation operators of the massless and massive quarks, and their Hermitian conjugate quantities are associated 
to the creation operators of the corresponding quarks. 

Considering these definitions H can be expressed in the traditional form 

^ = E E (^o(^) «^l^)«K^) + eoif) dlit)ciit))+ 
■f i=±i 

E E (^™/(^) bl{-^)hr{-^)+^ra,{-f) 4( j^)^^ ( jt)). (49) 

^ r=l,2 

The fields are now written as follows 

*(x) = ^„,(x) + *o(a:) = Y. *-(^)' (50) 



*o(^) = ^tE E {a\{-^)<{-^)^M-^ {e,{-f)t-^.l^))- 



^ t /=±i 



Ci(^)fo(^)exp(z (eo(y)t + ^.^))) 
Mo(a;) +wo(a;), (51) 

*„(a;) = TT^E E (^'■(^)<n,(^)exp(-i (e„,(^)t - ^.^)) + 

4(y)<n,(l^)exp(z (e™,(y)i + l^.^))) (52) 



^■^.=o 



Um{x) +Vm{x). 
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Now, in order to further define the quantization, it can be also taken into account that in the considered free theory, 
the expression for H should also be consistent with the assumed classical time dependence of the modes, which was 
employed to write the time derivatives of the field defining H . For the massless case, this condition imposes that 
the time evolution of the operators generated by the usual similarity transformation with the evolution operator of 
the operators ai{p) and cjip) should reproduce the classical time dependence of the field ^0(2;)- That is 

ai{p ,t) = exp(i H t) ai{p ) exp(— iiJ t). 

= exp(i eo(^) aj{~f)ai{~f) t) ai{~f) exp{-i eo{~f) a\{~^)ai{~f) t). 
= exp(-so i eo(lt) t) ai0) -)■ exp(i eoC^) t) aiCf) (53) 

c|("^,t) =exp(« H t) cjCj}) exp(-i H t). 

exp(i eo(y) c\{-f)ci{t) t) c\{-^) exp(-^ eoC^) c\{-^)ci{t)t)- 

= exp(i soeo(^) t) c\{~^) -^ exp(-i eoC^) t) q(^,0). (54) 

Thus, the Heisenberg time evolution of the operators will coincide with the classical one, upon the selection so = —1. 
This indicates that the massless field states should have negative norms. 

On the contrary, the similar consistency conditions directly lead to the usual selection s^ = 1 determining the usual 
positive norm for the one particle massive states. We thus conclude that, as suggested by the initial form of the 
free propagator, the Fock space of the free theory show an indefinite metric for fermion states, in which massless 
one particle states have negative norm and massive states show a positive one. 

C. The operator field equations 

Let us show that the defined field operator and its Hermitian conjugate solve the quantum version of the clas- 
sical Hamiltonian and Lagrange equations following by applying the standard quantization rules to the generalized 
Hamiltonian system. 

Since the entering two fields in any of the following expressions, are defined by two sets of creation or annihilation 
operators with a vanishing anticommutator between any element in one of the sets with any element pertaining to 
the other set, the following standard vanishing commutator relation at different spatial point and coincident time, 
follow for the massless and massive components of the fields 

{*o(a;),*o(a;')}Uo=xi =0, 
{*o(a;),*m(x')}Uo=4 =0, 
{vl/o(x),vI/t„(x')}U„=.i-0, 
{*™(x),*™(x')}Uo=.^=0. (55) 

Their Hermitian conjugate expressions furnish a similar set of relations including Hermitian conjugate fields. For 
the relation between the fields and their Hermitian conjugate (h.c.) being of the same massless or massive class, it 
follows 

{^I'o(x),vI/t(a:')}U„=,;,=ri+T2, 

exp {i {eoCf') - eo("^))xo - iCf' .^ - ~f .^)) 
= -TT^ui{p)uj{^) cxp { - i'f .{~^' -^)) 

{vI/o(a:), ^l{x')}U=., = - S^^H-^ ^V, (56) 

in which the term T2 is following after similar steps. The symbol ® will indicate the external product of two spinors 
leading to a spinor matrix. The negative metric anticommutation relation and spinor completeness relations had been 
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employed 



{4(^0, aK^)} = - ^M' S^^Kf'.t) I, (57) 

{4(^0, q(^)} = - Sij, S'-^'Hf^t) I, (58) 



i=±i 



,5(^)(^ - ^') = I] -^ exp ( - I ^.(^' - it)) 



The completeness condition follows from the massless spinor definitions in ([24|) and ([26| . Similarly, for the 
commutation relations between the massive field and its h.c. it is possible to write 

{^rn{x)M.i^')}\.o=.i,=Sl+S2, 

^1 = -^ E E Mp)®nUf'){bli-f'),Ki-f)}x 

exp(i {emfijt') - ejnf{jt))xo - i Cf'.lt' - jt.Ht)) 
= — ^ Ur{p)<^ul{']t)exp{-i^.{~^' ~1^)), 

52 = - X Vr{p)(E)vU~f)exp{-i^.{-^' ^it)), 

{*„,(x),vI.tja;')}Uo=.i ='5(^Hl^-^') ^- (60) 

In this massive case, the relations employed to arrive to the expression of the commutator were 

{bl,{-f'),K{-f)} = Sr,r'S'^''H'f',f) I, (61) 

{dl,if'),dr{f)} = Sr^r'S^^Hf'^lt) I, (62) 

r=±l 

The obtained commutation relations between the massless and massive components of the fields, then allow to write 
for the total fields ^ and 4*^, the commutations relations 

{^{x), *(x')}U„=x^ = {^\x),^H^')}U=.'o = 0' (64) 

{^{x),^Hx')}U,=,,„ - {*o(a;), *S(a;')}U„=.i + {^m{x), *L(x')}Uo=.i (65) 

= 0, 

which reproduce the similar relations between the Poisson brackets between the coordinates in the classical theory. 
Let us now represent the momenta operators defined in terms of the operators '^ (x) , 'i'' (x) by the same classical 
expressions 

n^{x) = -^¥{x)+Kdt^''{x)l3, (66) 

Il^f{x)^-'-^{x)-><Pdt^ix), (67) 

written in terms of the defined field operators. Then, for the anticommutator of them with the "coordinates" ^(x) 
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and ^(x) it follows 



= Ri + R2, 

Ri = -^ xJ2^o{~f){Mp) <» uj{~f) - vi{p) vl{~f)) exp ( - z -f.{-^' - ^)), 

i?2 = -^ ^H'^"(^)("'-(P) ® "'■(^) " ^'■(P) ® i'r('^))exp( - i f.i'^' - ^)). (68) 

But, after employing the spinors definitions p9)) . pOI) and p4)). (|26|) it follows 

/ ^\ 

^(z.z(p)®«|(l^)-t;,(p)®z;/(^))^U..^ 1^1 , (69) 

( I tjt \ 

r=i,2 em/IP J \ ^;r;- -J / 

which gives the commutation rule 

{vI/(:r),n*(a;')}Uo=xi = » '5^''^^ ^ ^V- (71) 

Therefore, the momentum operators defined in terms of the field operators as following from the classical equations 
also satisfies the usual quantization condition defining the anticommutator between the momenta and coordinates 
as given by the classical Poisson bracket between them times the imaginary unit. 

Considering that in identical form as it happens for the classical Grassman fields and momenta (n^(a;'))t = 
— n,j,t(a;'), and after taking the Hermitian transpose of ([71]) it also follows 

{'$^{x),U^,{x')}U,=,,^^^ S^''^-^ -^')I, (72) 

for the other set of conjugate pair of coordinate and momentum in the system. 

The anticommutation between momenta operators associated to the same fields ^I* or ^I^t 

{n*(a;),n*(x')}Uo=x;, ={n*t(a:),n^t(x')=0, (73) 

directly follow because any of the creation or annihilation operators included in the set S of them defining each of the 
entering momenta does not has its Hermitian transposed included in this same set S. 

Further, after employing the definitions of the momenta in terms of the fields, the anticommutator between momenta 
operators associated to the fields ^I' and 'i''' can be expressed in the form 

{n*(:r),n*t(x')}Uo=.;, = -x^dt^^ix') (3,(3 dt'${x)}+ 

'-^{{'$Hx'),(3dt^{x)}-{d,'i\x')(3,^{x)}) 

= -x^dt^l(x') (3,(3 dt^oix)} - x^dt^lix') P,(3dt^„,{x)}+ 

- '-{^Hx'),U^,ix)} ~ '-{U^ix'),^ix)}. (74) 

But, the direct evaluation leads to 

-^^{dt^lix')(3,(3dt^oix)} - I E ^ e^P(* ^-(^ - ^'))' (75) 



-K^dt^l{x')(3,(3dt^.m{x)} = -T7Y. 2 ^ exp(z ^.{1^ - ^0), 



V ^ mj 



(76) 



-'-{^Hx'),n^,{x)} - '-{n^{x'),^{x)} = ^ Ee^p(* ^-(^ - ^')), (77) 
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which after substituted show the resting standard anticommutation rules 

{nvp(x),n*t(x')}Uo=.;,=0. (78) 

Therefore, the new defined quantum fields ^(x) and ^^(x) and their momenta n^(x) and Il^t satisfy the standard 
quantization rules of the classical classical Grassman counterparts. The Hamiltonian of the system H can be now 
defined in terms of the classical one by substituting the classical fields and canonical momenta by the above defined 
operators in a way assuring that H is hermitian. The same classical expression gives a Hermitian representation in 
the form 

dx^{-{-n^,{x) /3n*t(x) + -*t(a;) f3^il{x)- 
X 4 

'-^\x) (i n*t(x) - ^n*(x) (i ■^{x))+ 

^*t(x)a' V,*(x) + X *t(x)^.^/3*(x)}. (79) 

Then, the use of the commutation relations 

{vI't(x),n^t(x')}U„=.^=*<5(^)(lJ^-li^')^ (80) 

{*(x), nv,(x')}U„=.^ = I s^°\i^ - 1^') I, (81) 

in conjunction with the rest of them among the fields and momenta which vanish, and the Heisenberg equation of 
motions for arbitrary operators A 

dtA{t)^t [H,A{t)], (82) 

permits to check that the fields and momenta obey the equations of motion 

dt^{x)^--P{n^,{x)+'-^{x)), (83) 

X 2 

dt^^^-{n^{x)+'-M'^{x))^, (84) 

K 2 

dtU^ix) =-^n^{x) p-^-^^x) p (85) 

2X 4:K 

+ i *tQ,«V, + >fV2*t^, (86) 

dtU^fix) = -^13 n*t(x) + —13 *(x) 

The further elimination of the momenta among these equations check that the "coordinate" fields also satisfy the 
original Lagrange equations 

{-i 7^a^ - m: o>2) *(x) = 0, (88) 

*+(x) («7^V^->< V2)=0. (89) 



V. SUMMMARY 

In this work, it has been studied the nature of the free quantization of the quark free Lagrangian associated to 
the recently proposed local and gauge invariant alternative for QCD for interacting massive quarks. It was argued 
that the free theory can be quantized in a way in which the massive quarks states can show a positive metric and 
the massless ones a negative norm. The quantization of the classical theory had been implemented, and an explicit 
representation of the fields and momenta operators solving the Heisenberg equations of motions was constructed. It is 
underlined that the accepted absence of asymptotic states for quark and gluons in QCD indicates that the appearance 
of indefinite metric in the quark sector of the theory does not necessarily present a definite limitation of the the 
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investigated alternative for QCD to describe nature. This conclusion follows because, precisely, the physical validity 
of the Lee- Wick theories rests in the absence of stable asymptotic states associated with the negative metric states 
due to the interactions [l3|- Therefore, the recognized today absence of free quarks and gluons in QCD can imply the 
satisfaction of this condition, through the confinement effect. The discussion also allowed to define a representation for 
the general classical action, and of the associated Dirac propagator of the Feynman expansion, which could improve 
the previous ones employed in Ref. [l|, for applications to high energy processes, where the massive quark states show 
approximate existence due to asymptotic freedom. It also could be expected to be of help in Nuclear QCD studies 
were also massive quark are relevant. 
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Appendix A 

In this appendix we sketch the derivation of the terms which define the standard formula for the energy operator 
as the superposition of the number operators for the various particles after multiplied by their respective energies. 
Below, the main algebraic steps in evaluating the general components of the energy operator (|37p are explicitly 
illustrated for two components. The evaluation for rest of the terms follows very similar algebraical lines. 

For the contribution of massless positive energy states it follows 



f dx'uiix)hid)uoix) = J2J2J2 E I 



X 

2 

■f, -f' i=±ii'=±i 



/3t''(^') r/3t''(y')) X 



^ jdx^eM^ {-^^^').-^)a\,{-t')ai{-t) 



■f, -f' l=±ll'=±l 

^.^)(1 - n') + l{l + 011^1)4 (P')a/(P) 
= EE^2/2|l^|at(lt)a,(l^) 

= EE ^oit)ali-f)M-f). (Al) 

-f i=±i 

In a very similar way for the negative energy contributions 

fdx\l{x)h{d)vo{x)=Y,Y. ^o{f)4itMt) (A2) 

•^ If i=±i 
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In the case of the terms related with massive quarks, the derivation of the energy term is sketched as foUows 



fdx'uiix) hid) u^{x) ^J2J2J2 J2 



emf{~f') + mf emf{~f)+mf 
z , — ^.. X -1 / : ;^ 



V 



^_ -fr r=l,2r' = l,2 V mf\i^ j y rnf\y 

^Jdx'^cMH'f-l^')-^) Xbl,{-f')brilt)x 

-^.y -^{eo{lt')eo{t}+f'-t} 



w 



\t) 



^r^-^(^)-'^(^) 



\ e™y(p)+m/ 

After evaluating the spinor products, expressing the spatial integral as the Kronecker delta of the two spatial 
momenta, and using the orthonormality of the two component spinors, this expression can be simplified in the way 
described below 



I dx^ulix) h{d) u^{x) ^J2J2 (^(^"V (^)^"/ (^) + ^•^)(1 - 7 ?^X ^ ^2 ' 

^■mf{p)+mf 2emf(P) 

+ m^ — j}'^ + 2m f eruf (1^) + mj 

{e„,,m+mfY )" 



^, '■=1,2 

^■mf{P)+mf 2e,nf{p) 

= J2e^^{t)emAt)bl{t) Kit)^^, 

-^ ^mf\P ) 

P, 

= E ^™/(p)&^(P )&'■(?) (A3) 

"^, r=l,2 



Finally, for the negative energy massive term, a closely similar calculation gives 



dx^vl^{x)h{d)v„,^\x)^ E e„,A~^) dli^) dr{-f). (A4) 

y, r-=l,2 
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